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Abstract 

We analyze the influence of the magnetic field generated by the supercurrents (self-field) on the 
current density distribution by numerical simulations. The thickness of the superconducting film 
determines the self-field and consequently the critical current density at zero applied field. We 
find an equation, which derives the thickness dependence of the critical current density from its 
dependence on the magnetic induction. Solutions of the equation reproduce numerical simulations 
to great accuracy, thus enabling a quantification of the dependence of the self-field critical current 
density with increasing film thickness. This result is technologically relevant for the development 
of coated conductors with thicker superconducting layers. 
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A decrease of the critical current density J c with increasing film thickness d was observed 
in thin superconducting films of YB^CusOy-^ on both single^ and polycrystalline^ - - sub- 
strates. Possible explanations are a degradation of the film or a change in the defect structure 
with thickness, but also the self-field of the sample causes a significant reduction if J c de- 
pends on the low magnetic induction B generated by the transport current^. Identifying the 
true cause of the decrease is vital for the development of coated conductors, which must be 
grown to higher thicknesses to enhance their performance. 

We employ numerical calculations of the current density distribution in a thin film to 
derive a practical approximation allowing the evaluation of the self-field depression of J c for 
a given J C (B) dependence. Therefore, it is essential to distinguish between the intrinsic J C (B) 
of the material, i.e., the (local) dependence of the critical current density on the magnetic 
induction B in the film, and the average critical current density as a function of the external 
applied field J C (H), which is determined in a transport measurement of the critical current. 
Note, that J C {B) provides information on pinning in the material and bridges theory and 
experiment. 

The procedure we employ is similar to the iterative algorithm used by Rostila et al.— At 
the beginning the external applied field determines the starting current density distribution 
J c (B(x,y) = HqH), which is constant across the cross-section of the film. After calculat- 
ing the self-field distribution B si (x,y) in the sample the current densities are updated to 
J c (B(x,y) = fioH + B s f(x,y)), which results in a new self-field distribution. This step is 
iterated until current density and magnetic induction satisfy the (arbitrary) material law 
J C (B) at every position in the film. For comparison with experiment J C {H) is computed by 
calculating the average current density flowing through the conductor at a range of external 
applied fields. 

Earlier 7 we used this procedure to extrapolate the power-law J C (B) oc B~ a dependence, 
which is commonly observed in thin films of YBCO at high fields, to the lowest fields and 
showed that despite the divergence of J C (B) the film carries a finite critical current at zero 
applied field because of the self-field of the sample. This result (see Fig. CP is particularly 
suited for our purpose, because the clear deviation between J C (B) and J C (H) at low fields 
emphasizes important features of the J C (H) curve: Starting from zero applied field, J C {H) 
remains approximately constant at the self- field critical current density J s f = J C (H = 0) 
as long as the applied field is negligible compared to the self-field of the sample. At high 
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FIG. 1. Approximations made to connect J s f(<i) and J C (B). The simulations of J C (H) (solid line) 
are carried out assuming J C {B) oc B~ 5 (dotted line) and result in a finite J s f at zero applied field. 
The simulated J C {H) remains approximately at this value up to applied fields of Bf , which is found 
by intersecting J s f (broken line) and J C (B), and approaches J C {B) above this field. 

applied fields, on the other hand, the transport current alters the field distribution only 
marginally and J C (H) is identical to J C (B). The applied field at which J C (H) becomes field 
dependent and approaches J C (B), can be estimated by comparison to the self-field of the 
film. We use iff = Bf/fiQ = 'yJgfd/Tc, which is equivalent to the field scale of thin films^, but 
we replace the constant J c of the Bean model by the average critical current density at zero 
applied field and include a factor 7, which is approximately 1.2 as determined graphically 
from Fig. [I] or by least squares fitting (see below). 

Because the sheet current density J s {d controls the self-field of the sample, the thickness 
dependence J s f(d) can be related to J C (B) by making the following approximations. We 
neglect any field dependence of J C (H) up to applied fields H = Hf, which leads to 

J C {H = 0) = J C (H = H { ) , (1) 
and assume that J C (H) and J C (B) are identical at applied fields above Hf. 

J C (H = H { ) = J C (B = B { ) . (2) 
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FIG. 2. Comparison of the simulated J s {(d) depression (points) to the solution of Eqn. 0] (broken 
lines) assuming a power-law J C (B). After fitting 7 ~ 1.2 the solution of the implicit equation is in 
excellent agreement with the numerical computation. (The a = 0.5 data is shifted up by a factor 
of 2 for clarity.) 

Combining both equations 



J C (H = 0) = J C (B = B t ) 
and inserting the definition of Bf results in the implicit equation 
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J sf = J C (B = /i 7^/vr) , 7 ~ 1-2 ■ (4) 

A graphical representation of the derivation is depicted in Fig. [1] and 7 is determined 
by intersecting J s f with J C (B). The final result Eqn. H] allows us to calculate the self-field 
critical current density as a function of thickness J s f(d) for any given J C (B). 

We test Eqn. H]by comparing its solution for a certain J C (B) to simulations of J s i{d) on a 
100 \xm wide film having a thickness between 100 nm and 3 (xm (the typical thickness range 
in experiments). Inserting, for example, the power- law 

UB) = h (-§-) ° , (5) 
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where J\ is the current density at B = Bi, generates a power-law also in the thickness 
dependence: 



Figure |2] compares simulations and Eqn. El We find by least squares fitting 7 R3 1.2, 
in agreement with the above result. Note further, that this pre-factor represents only a 
constant vertical shift in Fig. [2] and that the dependence on the thickness, which is in 
excellent agreement with the numerical simulations, is thus entirely the result of Eqn. HI 

After confirming the analytical solution for a power-law J C (B) we can analyze experi- 
mental data. Fits to J s {(d) curves are, for example, available from Ijaduaola et al. 5 , who 
finds J s {{d) oc d~ 0A in all but the thinnest film. According to Eqns. |5] and [6] the exponent 
a of J C (B) translates into an exponent of a/(l + a) in the thickness dependence. We infer 
a = 0.6, which is reasonably close to a ~ 5/8 determined in measurements of J C (H) at 
applied fields much above the self-field, where J C (B) and J C {H) are identical. The J s {{d) de- 
pression observed in this work can therefore be explained without any additional assumption 
by continuing the power-law J C (B) down to the low self- fields of the sample. 

Another popular parametrization of J C {B) is a generalized form of the Kim model 



which reproduces a power-law if B ^> Bq and is constant if B <C Bq. The maximum cur- 
rent density of J C (B = 0) = J limits the self-field of the sample to about 1.2fi J Q d/7T, which 
is roughly 14 mT at 3 (xm (the thickest film of the simulations) if we assume Jo = 10 10 Am" 2 . 
Thus we expect that the self-field, which increases linearly with thickness, significantly af- 
fects J s f{d) only if it is comparable to or exceeds Bq. Calculations for three different values 
of Bq are given in Fig. [3] and show the expected behavior: The self-field of the thinnest 
sample is approximately half of Bq = 1 mT and causes a noticeable J c decrease, whereas the 
self- field of the sample is insignificant and J c (d) approximately independent of thickness, 
if Bq = 100 mT. The intermediate B = 10 mT shows a transition from the first to the 
second case; in the thinnest films J s {{d) is nearly constant and the depression evolves as the 
thickness (and the self-field) increases. The (numerical) solution of Eqn. H] accounts for all 
three situations and matches perfectly the simulations using again 7 w 1.2. 
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FIG. 3. Same as in Fig. [2j but for a generalized Kim-model. Depending on the relative strength of 
Bq to the self- field of the sample, J s f(d) is almost constant if Bq B{ (squares) or decays rapidly 
if Bq <C Bf (circles). Equation |4] fully accounts for this behavior. 

We demonstrated how to calculate J B f(d) if a theoretical J C (B) is available, but our re- 
sult is also technologically relevant for investigating the reason of the J s f decrease in coated 
conductors with thicker superconducting layers. For such an analysis J C (H) should be mea- 
sured on the thinnest sample, where substrate interface effects do not play a role anymore. 
From this data J C {B) is determined down to inductions of about Bf of the thinnest film. 
The remaining uncertainty in J C (B) (we suggest to assume a constant J C (B) at inductions 
below Bf if no theoretical model is available) will introduce only a minor error in calculating 
J s f{d) for thicker films, provided they carry higher sheet current densities and thus elevate 
Bf above the value of the thinnest film. 

We wish to emphasize, that the effect of the self-field is not an alternative explanation 
for a thickness dependent critical current density and will always be present, if the critical 
current density depends on the magnetic induction. A depression of J s {{d) significantly 
different from that calculated from J C (B) is thus indicative of an additional mechanism, 
such as a change in the material or pinning properties with thickness, determining the 
performance of the thicker conductor. If, however, the calculated and the measured J s f(d) 
are similar, the self-field controls the thickness dependence. In this case an improvement in 
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pinning in the entire sample is necessary to enhance J C {B) and to improve the conductor. 

In summary we have derived a practical equation, which bridges theory and experiment 
by relating the (local) dependence of the critical current density on the magnetic induction 
J C {B) to the thickness dependence of the average critical current density at zero applied field 
J s f(d). Solutions of the equation correctly predict the thickness dependence for two common 
J C {B) models and excellent agreement is reached by fitting a single pre-factor. Thus the 
influence of J C (B) on the measured J B {(d) in transport experiments can be quantified, which 
is technologically important, because it allows investigating the depression of the critical 
current density observed in coated conductors when they are grown to higher thicknesses. 
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